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SUMMARY 

Dynamic  fracture  and  crack-arrest  responses  of  the  well-discussed  wedge- 
loaded  rectangular,  double  cantilever  beam  (WL-RDCB)  specimens  machined  from  a 
3.2  mm  thick  annealed  polycarbonate  sheet  were  studied  by  dynamic  photoelasticity 
and  dynamic  finite  element  analysis.  A  small,  but  finite  plastic  yield  zone 
surrounded  the  propagating  crack  which  resulted  in  crack  velocities,  dynamic 
fracture  toughness  and  crack  arrest  toughness  lower  than  those  observed  previously 
in  6.4  mm  thick  polycarbonate  fracture  specimens. 


INTRODUCTION 

Current  experimental  and  numerical  investigations  in  dynamic  fracture  involve 

extensive  use  of  dynamic  photoelasticity  [1,2]  dynamic  caustics  [3],  dynamic  finite 

difference  analysis  [4,5],  and  dynamic  finite  element  analysis  [6,7].  The  results 

published  in  the  above  are  based  on  the  theory  of  elasto-dynamlcs  and  have  evoked 

a  controversy  regarding  the  existence  of  a  material  property,  i.e.,  a  unique 

» 

dynamic  fracture  toughness,  K^,  versus  crack  velocity  relation,  a,  as  a  material 
property  which  governs  dynamic  crack  propagation  and  arrest  [1,2,6,8].  Despite 
this  unresolved  controversy  on  elastodynamic  fracture,  the  increasing  research 
efforts  in  this  field  should  eventually  lead  to  a  phenomenological  dynamic  frac¬ 
ture  criterion  which  can  be  used  in  either  finite  difference  or  finite  element 
codes  to  predict  crack  propagation  histories  in  actual  structures.  For  example, 
recent  numerical  simulation  of  dynamic  fracture  in  a  thermally  shocked  structure 
[9]  correlated  well  with  experimental  results  [10]  and  undoubtedly  such  practical 
applications  will  Increase  as  the  controversies  in  elastodynamic  fracture  are 
resolved. 

Unlike  the  efforts  in  elastodynamic  fracture,  little  basic  research  in  dynamic 
ductile  fracture,  particularly  in  the  presence  of  large  scale  yielding  has  been 
conducted  to  date.  Earlier  Investigations  on  dynamic  ductile  fracture,  including 
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those  of  Kanninen  et  al  [11],  Hahn  et  al  [12]  and  Ogasawara  et  al  [13]  as  well 
as  the  recent  papers  by  Koshiga  et  al  [14]  and  Kora  et  al  [15]  involved  pheno¬ 
menological  modelings  of  specific  experimental  observations.  While  these  results 
are  useful  in  predicting  dynamic  ductile  fracture  responses  in  specific  structures, 
they  fail  to  address  the  basic  laws  governing  dynamic  ductile  fracture. 

One  major  difficulty  in  dynamic  ductile  fracture  research  is  the  lack  of  suit¬ 
able  experimental  technique  with  which  the  dynamic  states  of  stress,  strains, 
plastic  yield  zone,  etc.  can  be  measured  directly  for  detailed  scrutiny.  In 
contrast,  dynamic  photoelasticity  and  dynamic  caustics  provide  optical  imaging 
from  which  the  dynamic  stress  intensity,  K^n  ,  among  others  can  be  extracted 
for  elastodynamic  fracture  studies.  Lacking  a  pure  experimental  approach,  one 
method  which  can  be  used  is  the  combined  experimental-numerical  procedure  for 
computing  dynamic  fracture  parameters.  In  this  procedure,  experimental  results 
such  as  crack  position  versus  time  in  a  fracturing  specimen,  is  used  to  drive  a 
crack  in  a  finite  element  model  of  the  fracture  specimen.  The  numerical  results,, 
such  as  the  equivalent  dynamic  stress  intensity  factor  for  elastic  analysis,  dyna¬ 
mic  J-integral  and  crack  opening  angle,  COA,  among  others  can  then  be  studied 

as  candidate  dynamic  ductile  fracture  criteria.  The  soundness  of  this  procedure 
when  used  in  elastic  analysis  has  been  verified  by  one  of  the  authors  [2,16].  It 
was  then  used  in  a  series  of  sensitivity  studies  on  the  influence  of  various 
dynamic  fracture  toughness,  KIQ,  versus  crack  velocity,  a,  relations  [17,18].  An 
obvious  extension  of  this  elastic  procedure  is  its  use  in  elasto-plastic  analysis. 

This  paper  reports  some  preliminary  results  obtained  in  the  course  of  deve¬ 
loping  an  experimental-numerical  procedure  for  determining  dynamic  ductile  frac¬ 
ture  toughness  which  could  be  used  to  characterize  the  dynamic  ductile  fracture. 

EXPERIMENTAL  ANALYSIS 

The  concurrent  experimental  analysis  in  this  investigation  consisted  of 
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dynamic  photoelasticity  using  thin,  i.e.,  3.2  mm  thick,  polycarbonate  wedge- 
loaded,  rectangular  double  cantilever  beam  (WL-RCB)  specimens  which  la  shown 
by  Figure  1.  The  annealed  polycarbonate,  from  which  the  specimens  were  machined, 
has  a  well-defined  yield  point  of  44  MPa,  flows  under  high  stresses,  and  exhibits 
tensile  instability  with  accompanying  Lueder 's  bands.  Under  high  strain  rate 
loading,  such  as  in  the  vicinity  of  a  propagating  crack  tip,  the  strain-rate 
sensitive  polycarbonate  fractures  in  a  cleavage  mode  without  the  large  plastic 
zone  ahead  of  the  propagating  crack  tip.  Such  ductile  to  brittle  transition  In 
dynamic  loading  of  polycarbonate  is  discussed  in  Reference  [19].  Despite  this 
cleavage  fracture  appearance,  the  dynamic  photoelastic -patterns  in  Figure  2(a) 
and  2(b)  show  that  the  propagating  crack  tip  is  preceded  by  a  small  yield  zone 
of  1.2  mm  and  0.4  mm  diameters,  respectively,  and  thus,  represent  dynamic  frac¬ 
ture  in  the  presence  of  small  scale  yielding.  The  elastic  region  surrounding  this 
small  yield  zone  also  provides  the  elasto-dynamic  state  which  can  be  used  to  ex¬ 
tract  apparent  elastic  fracture  parameter  for  characterizing  such  dynamic  ductile 
fracture.  For  example,  the  isochromatics  form  the  standard  loops  from  which  appa¬ 
rent  dynamic  stress  intensity  factor  can  be  extracted  through  the  use  of  well- 
knpwn  elasto-dynamic  near-field  solutions  [20].  Possible  errors  involved  in 

characterizing  the  elasto-plastic  dynamic  field  with  such  straight  forward  appli- 

* 

cation  of  elastic  analysis  will  be  discussed  later.  The  16-spark  gap  Cranz-Shardin 
camera  and  the  associated  dynamic  photoelastic  systems  used  to  record  16  discrete 
dynamic  photoelastic  patterns  in  the  fracturing  polycarbonate  WL-RDCB  specimen  has 
been  described  in  many  previous  papers  and  thus  will  not  be  repeated  here. 

In  the  as-received  condition,  the  polycarbonate  sheet  exhibits  considerable 
residual  stress  distribution  and  thus  the  sheets  were  annealed  after  rough  cutting. 
This  annealing, which  caused  some  distortion  and  shrinkage, consisted  of  overnight 
heating  at  160°C,  followed  by  gradual  cooling  at  the  rate  of  5*C  per  hour.  The 
starter  crack  consisted  of  0.4  mm  wide  edge  crack  approximately  25  mm  length  which 
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was  machined  and  then  chiseled  to  simulate  a  somewhat  blunt  crack  tip.  Static 
and  dynamic  stress-fringe  constant  and  Poisson's  ratio  at  various  strain  rates, 
all  of  which  were  determined  previously  [21],  were  used  in  data  reductions. 
Static  modulus  of  elasticity  of  E  -  2.00  GPa  and  the  yield  strength  of  44  MPa 
were  determined  separately  for  this  particular  polycarbonate  sheet. 


CRACK  TIP  ISOCHROMATICS  IN  THE  PRESENCE  OF  SMALL  SCALE  YIELDING 

The  influence  of  small  scale  yielding  on  the  dynamic  isochromatics  surround¬ 
ing  the  crack  tip  was  studied  by  a  static  Dugdale  strip  yield  model  (DM)  as  shown 
in  Figure  3.  Since  the  Uestergaard's  stress  function  for  a  partially  pressurized 

semi-infinite  crack  tip  is  known  [22],  the  crack  tip  Isochromatics,  t  ,  for  a 

in 

stationary  crack  with  a  Dugdale  strip  yield  of  length  ry,  can  be  expressed  as 

t  “-r[(o  -o  )^  +  4t^  ]2  (la) 

m  2  xx  yy  xyJ 


where 


o  -a  ■  Kx 
xx  yy  I 


7  i  T*  +  sin 

,ln~ 

C-)!*  1  *  co.e 

rDM  ™ 


2t  - 
xy 


r  0_„  30_m 

y  DM  ,  DM 

—h,  7~ 

>2  +  1  +  2  ^,C“e™ 


where  r^  and  0^  are  the  polar  coordinates  with  origin  at  the  DM  crack  tip  in 
Figure  3.  The  stress  intensity  factor.  Kg,  corresponds  to  the  negative  of  the 
stress  intensity  factor  at  the  DM  crack  tip  without  the  Dugdale  strip  yield  zone 
From  the  equality,  the  length  of  the  DM  zone  becomes 
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where  is  the  yield  strength  of  the  material. 

The  resultant  crack  tip  isochromatics  in  the  presence  of  a  DM  zone  can  be 
obtained  numerically  for  the  highly  loaded  crack  tip  at  the  onset  of  crack  pro¬ 
pagation.  For  estimated  stress  intensity  factor  of  ■  2.44  MPai4i,  yield  strength 

of  a  *  44  MPa  and  an  lsochromatic  f range  order  of  7-1/2  or  t  >84.8  MPa  in  the 
y  max 

3.2  mm  WL-RDCB  specimen.  The  DM  strip  yield  length  is  ry  -  1.3  mm.  Figure  4 
shows  the  estimated  crack  tip  isochromatics  with  the  DM  zone  and  the  experimentally 
observed  isochromatics  of  Figure  2(a).  The  significant  difference  between  the 
estimated  and  experimental  isochromatics  indicate  the  amount  of  plasticity  correc¬ 
tion  necessary  to  increase  the  elastic  Kj  in  the  presence  of  small  scale  yielding. 

The  procedure  then  for  estimating  the  apparent  stress  intensity  factor 
from  the  isochromatics  in  the  presence  of  small  scale  yielding  is  to  first  estimate 
the  location  of  the  far  end  of  the  crack  tip  yield  zone  or  the  DM  crack  tip  loca¬ 
tion  in  the  photoelastic  pattern.  An  overdetermlnistlc  least  square  method  [23]  is 
then  used  to  fit  at  the  DM  crack  tip  a  two  parameter,  elastic  crack  tip  stress 
field  to  the  experimentally  recorded  isochromatics.  Using  the  elastic  stress  inten¬ 
sity  factor  thus  computed,  the  DM  strip  yield  length  is  determined  by  equation 
(2)  and  the  estimated  lsochromatlcs  is  reconstructed  numerically  by  the  use  of 
equation  (1).  Should  the  computed  r^  be  substantially  different  from  the  measured 
r^,  then  r^  must  be  remeasured  and  the  above  numerical  procedure  is  repeated. 
Fortunately,  such  repeated  measurement  was  not  necessary  in  the  test  analyzed  for 
this  paper.  The  ratio  of  the  maximum  radial  distances,  rM  and  in  Figure  4,  is 
the  plasticity  correction  factor.  This  procedure  must  be  repeated  until  the  suc¬ 
ceeding  plasticity  correction  factor  approaches  unity.  However,  in  view  of  the 


approximate  nature  of  the  entire  procedure,  this  iteration  procedure  was  not 
programmed  into  this  data  reduction  scheme. 

The  above  procedure  is  based  on  static  crack  tip  yielding  and  obviously  must 
be  modified  for  dynamic  analysis.  It,  however,  provides  a  first  order  estimate 
of  the  plasticity  correction  involved  in  deducing  the  apparent  dynamic  stress 
intensity  factor  from  the  dynamic  isochromatics  in  the  presence  of  small  scale 
yielding. 

DYNAMIC  FINITE  ELEMENT  ANALYSIS 

The  dynamic  finite  element  algorithm  used  in  this  study  differs  with  that 
presented  two  years  ago  [24].  The  new  fracture  package  [2]  includes  a  mechanism 
for  gradual  release  in  crack  tip  nodal  force  as  well  as  an  iteration  routine  to 
assure  agreement  between  the  prescribed  and  computed  crack  tip  nodal  forces  in 
this  explicit  code.  Figure  1  shows  the  finite  element  breakdown  used  in  this  plane 
stress  analysis  which  is  considered  more  appropriate  for  these  thin  fracture  spe¬ 
cimens,  despite  the  cleavage  fracture  surfaces.  The  dynamic  finite  element  code 
was  executed  in  its  elastic-plastic  mode.  The  coarseness  of  finite  element  break¬ 
down,  i.e. ,  3.2  mm  square  in  comparison  to  the  observed  maximum  plastic  yield  rone 
of  1  mm,  however,  forced  the  actual  computation  into  its  elastic  mode  except  for 
the  initial  1-2  microseconds  in  crack  propagation  time.  As  a  result,  the  dynamic 
stress  Intensity  factor,  K^n  ,  was  computed  elastically  using  the  dissipated 
energy  at  the  released  node  with  Freund's  equation  [25]. 

RESULTS  OF  DYNAMIC  PHOTOELASTICITY 

A  total  of  seven  dynamic  photoelastic  tests  were  conducted.  Figure  2  shows 
two  typical  dynamic  photoelastic  patterns  of  a  fracturing  polycarbonate  WL-RDCB 
specimen.  These  dynamic  photoelastic  patterns  show  crack  tip  shadow  patterns  which 
are  believed  to  be  caustics  [3]  generated  by  the  large  plastic  strains  in  the  con- 
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fined  plastic  region  ahead  of  the  propagating  crack  tip.  The  size  of  such  a 
plastic  zone  ranges  from  ry  ■  1.3  n  to  0.04  mm.  Minor  distortions  of  the 
otherwise  elastic  isochromatics  in  the  vicinity  of  this  plastic  region  is  also 
noted. 

Apparent  dynamic  stress  intensity  factor,  K^yn  ,  was  computed  by  using  the 
data  reduction  procedure  described  previously.  Figure  5  shows  the  apparent 
dynamic  fracture  toughness,  which  were  obtained  by  this  procedure,  during  dynamic 
fracture  and  crack  arrest  in  a  polycarbonate  WL-DCB  specimen.  The  general 
variation  of  these  experimental  results  differs  slightly  with  those  obtained  by 
others  [3,7].  The  apparent  dynamic  fracture  toughness  continues  to  decrease  in 
the  former  while  distinct  regions  of  contant  KID  after  crack  propagation  were  ob¬ 
served  in  the  latter.  The  crack  initiation  stress  intensity  factor  of  - 
7.37  MPav^n  is  more  than  twice  of  the  pop-in  fracture  toughness  of  ■  3.4  MPai/m 
measured  previously  [21].  The  measured  apparent  crack  arrest  fracture  toughness 
of  KIa  »  1.0  MPa  m,  however,  is  considerably  lower  than  this  pop-in  fracture  tough 
ness  and  lower  than  the  KJa  »  1.4  MPaf'm  reported  for  the  6.4  mm  polycarbonate 
fracture  specimens  [2]. 

Also  shown  in  Figure  5  is  the  crack  velocity  change  with  crack  extension  in 
this  WL-RDCB  specimen.  This  crack  velocity  history  was  duplicated  in  the  other 
six  specimens  and  are  similar  in  shape  to  those  reported  by  others  [3,26]  but  the 
terminal  velocity  is  lower  than  that  in  Reference  [2]. 


RESULTS  OF  DYNAMIC  FINITE  ELH1ENT  ANALYSIS 

A 

Figure  5  shows  the  calculated  dynamic  stress  intensity  factor  obtained 
through  elasto-dynamlc  finite  element  analysis.  As  mentioned  previously,  the 
finite  element  analysis  was  executed  in  its  elastic-plastic  mode,  but  elastic 
analysis  prevailed  except  for  the  initial  1-2  microseconds  of  crack  propagation. 

I 

Agreements  between  calculated  dynamic  stress  intensity  factors  and  the  measured 
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dynamic  fracture  toughness  during  crack  propagation  history  are  apparent.  The 
precipitous  drop  in  calculated  dynamic  stress  intensity  factor  at  the  initial 
phase  of  crack  propagation  can  be  in  part  attributed  to  the  infinite  plate 
syndrome  discussed  by  Perl  [27].  It  was  also  found  that  this  precipitous  drop 
in  dynamic  stress  intensity  factor  could  be  moderated  by  assigning  a  gradually 
increasing  crack  velocity  during  the  first  30-50  microseconds  of  crack  propaga¬ 
tion  following  Reference  [27].  Recent  experimental  results  [20,28,29],  however, 
show  that  the  crack  velocity  of  a  blunt  starter  crack  immediately  reaches  the 
terminal  velocity  upon  propagation  and  does  not  gradually  approach  the  terminal 
velosity  as  postulated  in  Reference  [27]  during  this  crucial  initial  period  of 
crack  propagation.  These  experimental  evidences  thus  precluded  adjustments  of 
crack  velocity  to  reduce  the  precipitous  drop  in  R^n  . 

Figure  6  shows  the  energy  partitions  during  crack  propagation  in  the 
WL-RDCB  specimens  discussed  previously.  Although  the  dissipated  plastic  energy 
at  the  plastic  yield  zone  was  not  confuted  in  this  coarse  finite  element  model, 
the  decreasing  total  energy  in  Figure  6  indicated  qualitatively,  the  small 
plastic  energy  dissipated  in  the  fracture  process  of  polycarbonate  WL-RDCB  specimen. 


CONCLUSIONS 

1.  Small-scale  yielding  at  the  propagating  crack  tip  does  change  the  shape  of 
the  crack-tip  dynamic  lsochromatlcs. 

2.  Crack  velocity,  dynamic  fracture  toughness  and  crack  arrest  toughness  are 
lower  in  the  presence  of  small-scale  yielding  at  the  propagating  crack  tip. 

3.  Dynamic  stress  intensity  factor  obtained  by  elasto-dynamlc  finite  element 
analysis  are  in  good  agreement  with  the  apparent  dynamic  fracture  toughness 
obtained  through  an  elasto-plastlc  data  reduction  of  the  dynamic  lsochromatlcs 
in  the  presence  of  small  scale  yielding.  Elasto-plastlc  dynamic  finite  element 
analysis,  does  not  appear  necessary  for  such  dynamic  fracture  problems . 
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4.  The  decreesing  total  energy  indicate  the  influence  of  diealpated  plastic 
energy  in  the  preaence  of  a  small  plastic  tone. 
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